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ABSTRACT 

In this paper we describe and evaluate new spectral line polarisation observations ob- 
tained with the goal of mapping the surfaces of magnetic Ap stars in great detail. One 
hundred complete or partial Stokes IQUV sequences, corresponding to 297 individual 
polarised spectra, have been obtained for 7 bright Ap stars using the ESPaDOnS and 
NARVAL high resolution spectropolarimeters. The targets span a range of mass from 
approximately 1.8 to 3.4 M , a range of rotation period from 2.56 to 6.80 days, and a 
range of maximum longitudinal magnetic field strength from 0.3 to over 4 kG. For 3 
of the 7 stars, we have obtained dense phase coverage sampling the entire rotational 
cycle. These datasets are suitable for immediate magnetic and chemical abundance 
surface mapping using Magnetic Doppler Imaging (MDI). For the 4 remaining stars, 
partial phase coverage has been obtained, and additional observations will be required 
in order to map the surfaces of these stars. The median signal-to-noise ratio of the 
reduced observations is over 700 per 1.8 kms -1 pixel. Spectra of all stars show Stokes 
V Zeeman signatures in essentially all individual lines, and most stars show clear 
Stokes QU signatures in many individual spectral lines. The observations provide a 
vastly improved data set compared to previous generations of observations in terms 
of signal-to-noise ratio, resolving power and measurement uncertainties. Measurement 
of the longitudinal magnetic field demonstrates that the data are internally consistent 
within computed uncertainties typically at the 50 to lOOer level. Data are also shown 
to be in excellent agreement with published observations and in qualitative agreement 
with the predictions of published surface structure models. In addition to providing 
the foundation for the next generation of surface maps of Ap stars, this study estab- 
lishes the performance and stability of the ESPaDOnS and NARVA1 high-resolution 
spectropolarimeters during the period 2006-2010. 
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* Based on observations obtained at the Canada-France-Hawaii 
Telescope (CFHT) which is operated by the National Research 
Council of Canada, the Institut National des Sciences de l'Univcrs 
of the Centre National de la Recherche Scientifique of France, and 
the University of Hawaii. Also based on observations obtained 
at the Bernard Lyot Telescope (TBL, Pic du Midi, France) of 
the Midi-Pyrenees Observatory, which is operated by the Institut 
National des Sciences de l'Univers of the Centre National de la 
Recherche Scientifique of France. 



1 INTRODUCTION 

The classification Ap identifies a (main sequence) A or B 
type star which displays peculiar chemical abundances, usu- 
ally combined with an observable magnetic field. Although 
other classes of chemically peculiar stars exist (e.g. Am stars, 
Hg-Mn stars, He- weak stars), these stars have been demon- 
strated to lack strong, organised magnetic fields at their sur- 
faces (e.g. Shorlin et al. 2002, Wade et al 2006, Makaganiuk 
et al. 2011). Ap stars appear to be the only class of middle 
main-sequence stars for which, in all cases, an observable 
magnetic field is present (Aurierie et al. 2007). 
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Since their discovery by Babcock in 1947, the magnetic 
fields of Ap stars have been established through observa- 
tion to have important global dipole components with po- 
lar strengths ranging from hundreds to tens of thousands of 
gauss. The symmetry axis of the dipole component is almost 
always significantly tilted relative the stellar rotation axis. 
In addition, Ap stars generally spin much more slowly than 
non-peculiar stars of similar masses (Stepien 2000), and as 
they spin they exhibit line profile variations attributed to 
rotational modulation of patchy, non-axisymmetric lateral 
and vertical distributions of chemical abundance in their 
photospheres. The distributions of abundance vary signifi- 
cantly from element to element: some are distributed rela- 
tively uniformly, while others show strong contrast; some are 
distributed in relatively simple patterns, while others show 
complex distributions. While it is generally accepted that 
the fundamental mechanism responsible for the chemical pe- 
culiarities is microscopic chemical diffusion (as described by 
Michaud 1970), the origin of chemical patchiness, and the 
relationship to the magnetic field, is poorly understood. 

The earliest studies of the magnetic field geometries of 
Ap stars interpreted the rotational variations of their lon- 
gitudinal magnetic fields in the context of Stibbs' Oblique 
Rotator Model assuming a simple magnetic dipole field (e.g. 
Babcock 1947, 1951; Stibbs 1950). However, with the ac- 
quisition of increasingly sophisticated diagnostic data (the 
mean surface field (or mean field modulus), and high- 
resolution line profiles), it became clear that the large-scale 
field topologies exhibited important departures from the 
simple dipolar model (e.g. Preston & Sturch 1967, Preston 
1969, 1970, Landstreet 1970, 1988, 1989). 

Leroy and collaborators (Landolfi et al. 1993, Leroy et 
al. 1993, 1994, 1995ab, 1996, Bagnulo et al. 1995, Wade et 
al. 1996) systematically studied Ap stars using broadband 
linear polarisation measurements and models, constraining 
the transverse component of the magnetic field. Importantly, 
they found that differences between the observed linear po- 
larisation variations and those predicted by the simple dipole 
model could not be fully explained by abundance inhomo- 
geneities alone (Leroy et al. 1996). With this work, they 
established a modified dipolar model with a trend toward 
an outward expansion of the field lines over some parts of 
the magnetic equator, and showed the potential of linear 
polarisation for diagnosing small-scale structure of the mag- 
netic fields of Ap stars. Thus, the observations and mod- 
eling undertaken during the latter half of the 20th century 
allowed progress from a simple view of the magnetic fields 
of Ap stars to a relatively sophisticated picture in which 
fields were known to show both global-scale and local-scale 
departures from a simple dipole. 

Leroy et al. (1996) commented that high-resolution 
spectropolarimetry represented the next step in further- 
ing the study of the magnetic field geometry of Ap stars. 
Four years later, Wade et al. (2000a) published the first 
compendium of phase- resolved high-resolution spectropo- 
larimetric observations of Ap stars in both circular and 
linear polarisation. Using the MuSiCoS spectropolarimeter, 
R = 35, 000 Stokes IQUV spectra with a median S/N of 
300 (per 4.6 kms~ pixel) were obtained for 14 Ap stars. 
While the quality of the spectra was sufficiently good to 
show the shape and phase variation of all Stokes parame- 
ters in mean Least-Squares Deconvolved (LSD) line profiles, 



measurement in individual spectral lines was restricted to a 
few particularly strong lines, principally those of Fe II mul- 
tiplet 42. Nevertheless, the Stokes profiles of 53 Cam were 
used by Bagnulo et al. (2001) and Kochukhov et al. (2004) to 
evaluate published magnetic models developed based on less 
sophisticated data. Those authors found that models based 
on so-called "magnetic observables" (e.g Bagnulo 2000) led 
to derivation of surface magnetic field characteristics that 
were not consistent with the detailed Stokes profiles, and 
that both circular and especially linear polarisation profiles 
were required for realistic reconstruction of the field. 

Following these conclusions, Kochukhov et al. (2002) 
for a 2 CVn and Lueftinger et al. (2010) for HD 24712 em- 
ployed the new Magnetic Doppler Imaging technique (MDI), 
described by Piskunov & Kochukhov (2002) and Kochukhov 
& Piskunov (2002), to construct high resolution maps of the 
surface vector magnetic field maps using Stokes IV obser- 
vations and by preferring a global low-order multipolar field 
structure. Maps using Linear polarisation profiles (Stokes Q 
and U) made by Kochukhov et al. (2004) and Kochukhov & 
Wade (2010) for the Ap stars 53 Cam and a 2 CVn. These 
maps were distinguished from earlier models in that they 
were computed directly from the observed polarised line pro- 
files, making no a priori assumptions regarding the large- 
scale or small-scale topology of the field. The MDI surface 
magnetic field maps of both stars revealed that their mag- 
netic topologies depart significantly from low-order multi- 
poles. In particular, both studies concluded that while the 
global topology of the magnetic field was reasonably smooth, 
the strength of the field was quite patchy, indicating com- 
plex structure on relatively small scales. Simultaneous map- 
ping of the distributions of the surface chemical abundances 
of several elements was also performed, allowing a compari- 
son between the local field properties and local photospheric 
chemistry. 

It is important to note that the observational material 
used in the MDI studies of 53 Cam and a 2 CVn represented 
the best data sets obtained from several years of MuSiCoS 
observations. In those spectra the uniquely valuable Stokes 
Q and U Zeeman signatures were only clearly detectable 
in 3 strong lines, with a significance (i.e. amplitude divided 
by error bar) of 5 or less. The relatively low signal-to-noise 
ratio and resolving power achievable with the MuSiCoS in- 
strument led to some ambiguity in the field reconstruction, 
and limited the useful sample of stars to those with bright 
apparent magnitudes, strong fields and sharp lines. Because 
MDI exploits the indirect resolution of the stellar disc due 
to stellar rotation, this means that those stars best suited to 
reconstruction (relatively rapidly-rotating stars, with conse- 
quentially weaker Stokes profiles) were inaccessible to Mu- 
SiCoS. As a result, only an extremely limited range of stel- 
lar properties (rotation, mass, temperature, magnetic field, 
etc.) which may influence the phenomena of interest could 
be studied using the MuSiCoS data. 

To address outstanding questions surrounding the de- 
tailed magnetic structure of Ap stars and the effect of the 
magnetic field on atmospheric chemical transport processes, 
we have acquired new higher-resolution and signal-to-noise 
Stokes IQUV spectra of a small sample of well-studied mag- 
netic Ap stars using the new generation of high-resolution 
spectropolarimeters. In this paper we describe the observa- 
tions obtained. We demonstrate the stability of the instru- 
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OBSERVATIONS OBTAINED WITH 
ESPADONS AND NARVAL 
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Figure 1. Phase coverage obtained for the target stars. Symbols 
denote rotational phases at which observations were obtained. 



mentation during the 5 years of observation by evaluating 
the internal and external agreement of the data. We illus- 
trate the quality of the observed Stokes profiles, comparing 
with MuSiCoS results and demonstrating that they repre- 
sent a qualitative step forward in our ability to diagnose the 
magnetic structure of Ap stars. 



2 TARGETS 

Targets selected for this study are bright Ap stars demon- 
strated to exhibit strong Stokes profiles by Wade et al. 
(2000a). We attempted to select targets spanning a large 
range of stellar physical properties as well as field strengths 
and geometries. Target stars were generally required to have 
the following characteristics to be suitable for this study: 

• Well-determined rotation periods - This study requires 
phase-resolved timeseries observations. All targets must 
therefore have well-determined, unambiguous rotation pe- 
riods so that each phase is observed correctly. 

• Suitable projected rotational velocity - The projected 
rotational velocity must be neither too rapid (rapid rotators 
> 50 kms -1 typically have shallow Stokes profiles which 
are challenging to detect and interpret) nor too slow (< 2 
kms" 1 whilst MDI can be applied to stars with such small 
rotational velocities, the advantage of Doppler tomography 
is lost). 

• Strong magnetic fields and variability - Because MDI 
relies on both the shape and variation of line profiles to 
determine the geographic location of magnetic and chemi- 
cal features, target stars should display strong and variable 
Stokes profiles. 

Ultimately, 7 targets were selected for monitoring. The 
target list is shown in Table [1] 



Both the ESPaDOnS and NARVAL instruments consist 
of a bench mounted cross-dispersed echelle spectrograph, 
fibre-fed from a Cassegrain-mounted polarimeter unit. These 
instruments are designed to overcome the limitations en- 
countered with MuSiCoS, with improved resolution (R = 
A/AA ~ 65000), sensitivity (approximately 15-20% through- 
put) and wavelength coverage from 369-1048 nm (with gaps 
at 922.4 to 923.4 nm, 960.8 to 963.6 nm and 1002.6 to 1007.4 
nm). The ESPaDOnS (Echelle SpectroPolarimetric Device 
for the Observation of Stars) spectropolarimeter is installed 
at the 3.6 m Canada- France-Hawaii Telescope (CFHT), and 
the NARVAL spectropolarimeter is installed at the 2 m 
Bernard Lyot telescope at Pic du Midi observatory. ES- 
PaDOnS and NARVAL are essentially identical instruments, 
with NARVAL constructed based on the experience of ES- 
PaDOnS. 

The polarimetric unit (at the Cassegrain focus) allows 
two orthogonal states of a given polarisation (circular or lin- 
ear) to be are recorded throughout the entire spectral range. 
The polarimeters of ESPaDOnS and NARVAL are of a sim- 
ilar design to the "Sempol" visitor polarimeter (Donati et 
al. 2003) used on the AAT (Anglo- Australian Telescope). 
The polarimeter is split into two parts, the upper part is 
for guiding and calibration, containing the guiding camera 
(a commercial FLI MaxCam series CCD camera is used in 
NARVAL and as of semester 2011 A a QSI Imaging CCD 
camera is used in ESPaDOnS), an atmospheric dispersion 
corrector and calibration wheel. The lower part contains the 
Fresnel rhomb retarders which are used to perform the po- 
larimetric analysis. These optics are in 4 different drawers 
which contain the following (in order): a half- wave rhomb 
(consisting of a pair of quarter- wave rhombs), a quarter- 
wave rhomb, a second half-wave rhomb, and finally a Fabry- 
Perot wheel on one side and a Wollaston/wedge plate slide 
on the other side. To enhance the achromaticity of the phase 
delay, and to keep it at 90 ± 0.5° across the whole optical 
domain, the rhombs are coated with a thin layer of MgF2, 
yielding a performance significantly better than achromatic 
crystalline plates which can vary by about 20° from the re- 
quired quarter-wave retardance. Another advantage of Fres- 
nel rhombs is that they do not produce detectable spectral 
ripples due to the fact that the fringe spacing is on the order 
of the pixel size. The half-wave rhombs can rotate about the 
optical axis by a specified angle. The final component, the 
Wollaston prism, consisting of two orthogonal calcite prisms 
that are cemented together, acts as a polarising beamsplit- 
ter. The two beams of light from the beamsplitter are trans- 
mitted by some 30 m of optical fibre to the spectrograph. 
ESPaDOnS includes a fiber agitator, which shakes the opti- 
cal fiber to remove modal noise that may be present. 

The ESPaDOnS spectrograph unit consists of a dou- 
ble set of high-reflectance collimators cut from a single 680 
mm parabolic mirror, with a focal length of 1500 mm. The 
grating is a 79 gr/mm monolithic grating with a dimensions 
of 200 by 400 mm. The camera lens is a fully dioptric f/2 
388 mm focal length lens, with a 210 mm free diameter (7 
lenses in 4 blocks, one block being a 220 mm quadruplet). 
For cross-dispersing, a high dispersion prism made of a train 
of 2 identical PBL25Y prisms with an 35° apex and 220 mm 
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Table 1. Stars discussed in this paper, along with ancillary data: HD designation, other name, V magnitude, spectral type, projected 
rotational velocity, rotational period, radius, inclination, the number of longitudinal field measurements obtained, the median longitudinal 
field uncertainty obtained, the number of net linear polarisation measurements obtained. Projected rotational velocities (vsini) have 
been measured in this study, as discussed in the analysis section. The periods were obtained from various sources as described in the 
Sect 6. Masses are those reported by Kochukhov and Bagnulo (2006) and the average longitudinal magnetic fields are taken from the 
catalogue of stellar effective magnetic fields (Bychkov et al. 2003), using values obtained from Least-Squares Deconvolved profiles if 
available. Temperatures are those reported by Kochukhov and Bagnulo (2006) and stellar radii as reported in Leone et al. (2000) for HD 
32633, Pasinetti Fracassini et al. (2001) for HD 40312, Kochukhov and Wade (2010) for HD 112413 and Wade (1997) for the remaining 
stars. Inclinations are taken from Stepieh (1989) for HD 32633, Rice, Holmgren and Bohlender (2004) for HD 40312 and Kochukhov and 
Wade (2010) for HD 112413 and from Leroy et al. (1996) for the remaining stars. 
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cross section is used. Up until the 2011A semester the de- 
tector used with ESPaDOnS was a grade 1 EEV detector 
with 2K by 4.5K 0.0135 mm square pixels (known as EEV1 
at CFHT). This was replaced in 2011A with a new E2V 
detector (named Olapa at CFHT). 

The spectrograph unit is mounted on an optical bench, 
which is housed in a thermal enclosure found in the inner 
Coude room at CFHT. 

This configuration yields full spectral coverage of the 
optical domain (from grating order 61 centred at 372 nm 
to grating order 22 centred at 1029 nm) in a single expo- 
sure. In polarimetric mode this should, in principle, achieve 
a resolution in excess of 65,000, but due to a charge trans- 
fer efficiency issue with the EEV1 CCD detector, the true 
resolution varies from approximately 68,000 in the blue to 
61,000 in the red. The peak throughput of the spectrograph 
(with CCD detector) is about 40% to 45%, bringing the to- 
tal instrument peak efficiency to a level of about 15% to 
20%. 

The configuration of NARVAL is much the same with 
the exception of an 2.8 arcsec aperture pupil (versus 1.6 
arcsec for ESPaDOnS). NARVAL does however benefit from 
better spectrograph thermal stability (by approximately a 
factor of 10) than ESPaDOnS, due to the use of a double 
thermal layer enclosure. 

The resolving power and signal-to-noise ratio of both 
instruments vary with wavelength in a predictable manner, 
as illustrated in Fig. [5] which shows the signal-to-noise ra- 
tio as a function of wavelength for four observations of a 2 
CVn (two acquired with ESPaDOnS, and two with NAR- 
VAL). The variation in the spectral resolving power R as a 
function of spectral order for a selection of observing nights 
is shown in Fig. [3] This figure illustrates that the charac- 
teristics of the instruments do not vary significantly on the 
timescales relevant for this project (nights to years). Vari- 
ability in atmospheric conditions (e.g. seeing) may well be a 
dominant contributor to the scatter in resolving power. 

In terms of observations the following steps occur; the 
two output beams from the Wollaston prism, which have 
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Figure 2. Signal-to-noise ratio (per 1.8 kms 1 spectral pixel) as 
a function of wavelength for two observations of a 2 CVn obtained 
with ESPaDOnS (top curves, 02 March 07 (airmass = 1.104) and 
13 Jan 09 (airmass = 1.058) both total exposures of 120s) and 
two obtained with NARVAL (lower curves, 20 Dec 06 (airmass 
= 1.250) and 11 Jan 09 (airmass = 1.005) both total exposures 
of 240s). Solid vs dashed lines indicate observations taken on re- 
spective nights. 



been analysed into the two components of circular polarisa- 
tion (using the quarter-wave rhomb) or linear polarisation 
(using the half- wave rhomb, are then carried by the pair 
of optical fibres to the spectrograph where two interleaved 
spectra are formed. The / component of the stellar Stokes 
vector is formed by adding the two corresponding spectra, 
while the V, Q or U polarisation component is obtained 
from the ratio method as discussed by Bagnulo et al. (2009) . 
To minimise systematic errors due to small misalignments, 
differences in transmission, effects of seeing, etc., one com- 
plete observation of a star consists of four successive sub- 
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Figure 3. The mean spectral resolving power as a function of 
spectral order for both ESPaDOnS and NARVAL for a selection 
of observing nights (3 observations from each instrument spanning 
12 months). Error bars represent the standard deviation of the 
resolving power of the various observations within each spectral 
order. 



exposures; for the second and third, the waveplate settings 
are changed so as to exchange the positions of the two anal- 
ysed spectra on the CCD. This same procedure is used for 
all polarisation spectra, Stokes V, Q and U. 

Calibration of the instrument uses a combination of tho- 
rium/argon and thorium/neon lamps, with the lamp calibra- 
tions taken at the beginning and at the end of each night 
for a primary wavelength calibration. Telluric lines are then 
later used to perform a second wavelength calibration during 
the reduction process using libre-ESpRIT. Filters are used 
to minimise the blooming on the chip at the red end of the 
spectrum. Two tungsten lamps are utilised for the fiat fields 
frames, with one low intensity lamp being used with a red 
filter and the other lamp being higher intensity and used 
with a blue filter. 

The reduction of observations is carried out at the 
observatories using the dedicated software package Libre- 
ESpRIT, which yields both the I spectrum and the V cir- 
cular polarisation spectrum and/or QU linear polarisation 
spectra of each star observed. It is important to note that 
libre-ESpRIT automatically finds and removes continuum 
polarisation. In this work each reduced spectrum is nor- 
malised order-by-order using a FORTRAN code specifically 
optimised to fit the continuum of these stars. 

A diagnostic null spectrum called the N spectrum, com- 
puted by combining the four sub-exposures in such a way 
as to have real polarisation cancel out, is also calculated 
by Libre-ESpRIT (again, see Bagnulo et al. 2009 for the 
definition of the null spectrum). The N spectrum tests 
the system for spurious polarisation signals. The result of 
the reduction and normalisation procedure are continuum- 
normalised, one-dimensional spectra in the form of wave- 
length, I /I c , V/I c , Q/Ic, U/I c , two independent normalised 
A" spectra Ni/I c and A^/Zc, and an error bar (computed 
by propagating photon uncertainties through the reduction 
procedure), tabulated pixel- by-pixel. 

Weak signatures are visible in LSD A^ profiles associated 



with some Stokes Q and U spectra. Experiments conducted 
in order to diagnose and mitigate polarisation crosstalk (see 
Barrick et al., in preparation) indicate that these signatures 
are related to this phenomenon. However, data acquired 
during early NARVAL runs in 2006 exhibit substantially 
stronger signatures. These signatures are reported by TBL 
staff to result from problems with the coatings on the A/2 
rhombs in use in 2006. The rhombs were subsequently re- 
placed, and no similar strong signatures are detected in later 
observations. 

Based on our comparisons of Stokes Q and U profiles 
obtained with the two instruments during the course of our 
observing program, it is clear that these N signatures are not 
diagnostic of any detectable contamination of the associated 
Q and U profiles. 

In this study we have obtained 100 polarimetric se- 
quences corresponding to 297 individual polarised spectra. 
The observations were initially obtained in classical observ- 
ing mode, and later in service mode at both telescopes. In 
total 48 sequences were obtained with ESPaDOnS and 52 
with NARVAL. The log of observations is reported in Table 
[31 and the achieved phase coverage of the stellar targets is 
illustrated in Fig. [T] The resulting reduced spectra are illus- 
trated in Fig. [8] which demonstrates the quality of the data, 
with Stokes VQU signatures seen in many individual lines. 



4 CROSSTALK 

During the commissioning of ESPaDOnS in 2004 it was 
found that the instrument exhibited crosstalk between linear 
polarisation and circular polarisation (and vice versa). Due 
to the relatively strong circular polarisation in the spectral 
lines of our targets, contamination of the significantly weaker 
linearly polarised profiles is potentially a serious problem. 

The contributing optical component to this crosstalk 
identified initially was the collimating triplet within the po- 
larimeter unit. This component was replaced in June 2006, 
resulting in a reduction of the crosstalk to the 5% level 
(from an initial level of 10-15 %). After replacing the triplet 
lens again in October 2008, the crosstalk appeared to have 
been reduced to 2-3%, but still exhibited strong temporal 
changes. It was discovered in October 2009 that the atmo- 
spheric dispersion corrector (ADC) was an important and 
previously unrecognised source of crosstalk within the po- 
larimetric unit. The ADC was replaced in the fall of 2009, 
and since that time the crosstalk has been small and stable, 
with crosstalk from Stokes V into Stokes U at the 0.5% level, 
and no measurable crosstalk from Stokes V into Stokes Q 
(i.e. below ~ 0.1 — 0.2%). The procedure and results of the 
investigations of the crosstalk, are reported by Barrick et al. 
(in preparation). The evolution of the ESPaDOnS crosstalk 
with time, based on the results of Barrick et al., is illustrated 
in Fig. H 

Less extensive crosstalk monitoring has been performed 
with NARVAL. Results from September 2009 indicated that 
with the ADC in place, the crosstalk was 3.1% from Stokes 
V to Stokes Q, and below 0.1-0.2% from Stokes V to Stokes 
U. Without the ADC in place, the crosstalk was reduced to 
2.1% in Stokes Q, but increased to 1% in Stokes U (illus- 
trating that the NARVAL ADC also introduces crosstalk, 
but again that it is not the sole source). It is important to 
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note that at the current time, no tests of how this crosstalk 
changes with time have been made with NARVAL. 

Although the crosstalk is now below 1% in ESPaDOnS 
and probably around 2% in NARVAL, observations for this 
project have been acquired over the last 4 years and during 
some of this time the crosstalk may have been higher (and 
almost certainly was for ESPaDOnS). It is important to un- 
derstand how crosstalk could affect the Stokes Q and U sig- 
natures measured in spectral lines. A series of test was there- 
fore performed to evaluate the importance of this crosstalk. 

The first test was to examine the Stokes profiles of the 
Fe n A6149 line. This line has a relatively large Lande factor, 
but in the linear regime of the Zeeman effect it is predicted 
be purely circularly polarised as a consequence of the a and 
7r components of this line having the same strength and iden- 
tical splitting. In this case we would interpret any signal in 
Stokes Q and U as due to crosstalk from Stokes V. We have 
examined A6149 in our spectra, as well as spectra discussed 
by Barrick et al. (in preparation). In no case do we observe 
any significant signal in this line. 

Then, we carefully examined one of the crosstalk di- 
agnostic observations of the cool magnetic star 7 Equ ob- 
tained during CFHT engineering time. As described by 
Barrick et al. (in preparation), on-sky crosstalk diagnosis 
employs observations of slowly-rotating magnetic Ap stars 
observed in all Stokes parameters at two positions of the 
CFHT's Cassegrain bonnette to unambiguously measure the 
crosstalk into Stokes Q and Stokes U. The diagnostic obser- 
vation of 7 Equ (obtained in July 2009, with a peak S/N 
of over 1000) yielded relatively high crosstalk levels of 2.3% 
in Stokes Q and 5.1% in Stokes U. By comparing the ob- 
servations with and without crosstalk, it was found that the 
crosstalk contributions to Stokes Q is within the noise, whilst 
the crosstalk contributions to Stokes U appears to be slightly 
above the noise, as illustrated by Fig. [5] It is important to 
note that the removal of the crosstalk using techniques as 
described in Barrick et al. (in preparation) requires the ac- 
quisition of a second series of observations following a rota- 
tion of the Cassegrain bonnette, a procedure which is not 
standard practice during regular observations. 

7 Equ is an extremely sharp lined star, which means 
the crosstalk contribution will have a greater effect than 
it would in the broader lined stars studied in this paper. 
Nevertheless, the potential influence of the crosstalk on MDI 
will be investigated and discussed in a future paper (Silvester 
et al. in preparation). 

These comparisons demonstrate that the contribution 
of crosstalk to the Stokes profiles is below, or at most just 
above, the level of the noise of the best-quality observations 
of Ap stars acquired with these instruments. 
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Figure 4. The total crosstalk of Stokes V into Stokes Q and U 
of ESPaDonS as function of time, as reported by Barrick et al. 
(in preparation) . The dashed boxes indicate periods during which 
observations were obtained as part of this investigation. 
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Figure 5. A comparison between the measured crosstalk in 
Stokes Q and U and the null spectrum for observations of 7 Equ 
taken with ESPaDOnS (16th July 2009) (top row), with the ac- 
tual observations shown below (lower frame). In the strong, mag- 
netically sensitive line Fe II 5018 line, the crosstalk in Stokes Q 
is below the noise, where as the crosstalk in Stokes U is slightly 
above the noise. 



5 LONGITUDINAL MAGNETIC FIELD AND 
NET LINEAR POLARISATION 

To examine the self-consistency of the new polarisation spec- 
tra, as well as to evaluate their consistency with published 
magnetic data for our targets, we have measured the mean 
longitudinal magnetic field and net linear polarization us- 
ing Least-Squares Deconvolution (LSD). LSD (Donati et 
al. 1997 and Kochukhov et al. 2010) is a multiline anal- 
ysis method that produces mean Stokes I and V profiles 



using essentially all metallic lines in the stellar spectrum. 
It assumes that the observed spectrum can be represented 
as the convolution of a single mean line profile with an un- 
derlying spectrum of unbroadened metal and helium lines of 
appropriate wavelength, depth and Lande factor (the "line 
mask" computed using spectrum synthesis; e.g. Wade et al. 
2000a). 

The LSD model allows the computation of single, av- 
erage Stokes I and V line profiles, usually characterised by 
a signal-to-noise ratio significantly higher than that of indi- 
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Figure 6. Comparison between Stokes Q and U profiles of HD 32633 in the 5018 A line for ESPaDOnS/NARVAL on the left and 
MuSiCoS on the right. Rotation phase for each observation are indicated. 



vidual spectral lines, scaling roughly as the square root of 
the number of lines used. 

The projected rotational velocity v sini of each star was 
derived by fitting a selection of spectral lines (typically in 
the 4500 A region) with a synthesised spectrum created us- 
ing the Synth3 spectrum synthesis code (Kochukhov, 2007). 
Each line in the selection had vsini determined by using 
a \ 2 fitting function which is included as part of the Bin- 
mag spectral visualization tool. The resulting mean (and 
standard deviation) values one reported in Table [1] Line 
masks for this study were compiled using Vienna Atomic 
Line Database (VALD, Kupka et al. 1999) "extract stellar" 
requests, with effective temperatures (adopted based on the 
literature) of each target. It should be noted when creating 
each line mask, a constant log g of 4.0 was used for all stars. 
Input chemical abundances were determined by a rough 
abundance analysis, based typically on lines around the 4500 
A region. The abundance analysis was performed by using 
the Synth3 spectrum synthesis code (Kochukhov 2007) with 
ATLAS9 solar abundance model atmospheres used to create 
synthetic spectra that were compared directly to the obser- 
vations. The input abundances were then adjusted in the 
synthetic spectrum until a reasonable agreement was found 
with the observed spectrum. A comparison between of two 
synthetic spectra, one computed using solar abundances and 



one using the final determined abundances, is shown in Fig 
[7J for HD 32633. The final abundances were then adopted 
in the line mask creation. By doing so we ensure that any 
uncertainties caused by the mask when performing the LSD 
analysis are minimal. 

As discussed by Shorlin et al. (2002), the LSD S/N is 
only weakly sensitive to the line-depth cutoff employed to 
populate the mask. Following their results, we have chosen 
to employ a line-depth cutoff equal to 10% of the contin- 
uum. Imposing such a cutoff has a related advantage: be- 
cause weaker lines are less likely to have published experi- 
mental Lande factors (and to generally have more poorly- 
determined atomic data), we pre-filter our line list to (statis- 
tically) exclude those lines with the poorest data. In addition 
Balmer lines are removed from the mask and the mask is 
restricted to the ESPaDOnS/NARVAL spectral range. Ap- 
plication of LSD to the data yields a set of mean profiles 
(Stokes /, Stokes V and N) for each reduced spectrum. 

The mean longitudinal magnetic field Bi was computed 
from each LSD profile set. This quantity was evaluated by 
computing the first-order moment of the Stokes V profile in 
velocity according to: 



-2.14 x 10 



uj(v — vo)V(v)dv 
\gc J[l-I(v)]dv 



(1) 
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Stokes V Fit 


Null N 


Star 


Degree 


x 2 


x 2 


HD 4778 


2 


1.41 


2.68 


HD 32633 


3 


1.95 


5.42 


HD 40312 


2 


1.96 


0.51 


HD 62140 


3 


1.29 


4.78 


HD 71866 


4 


1.02 


3.70 


HD 112413 


3 


0.83 


1.06 


HD 118022 


1 


1.92 


1.12 



Table 2. Reduced x 2 s of harmonic fits to the longitudinal mag- 
netic field phase curves determined from NARVAL and ES- 
PaDOnS data 

(Mathys 1989, Donati et al. (1997), Wade et al. 2000a) where 
vo is the centre-of-gravity of the Stokes V profile, g is the 
integrated mean Lande factor and A is the weighted mean 
wavelength of all the lines included in the mask. LSD pro- 
files were locally re-normalised to a continuum level of 1.0 
before evaluation of Eq. (1). Uncertainties associated with 
B{ were computed by propagating the formal uncertainties 
of each LSD spectral pixel through Eq. (1). LSD profiles are 
extracted for each star uniformly weighted to the same lande 
factor, line depth and wavelength. 

To determine the net linear polarisation (see e.g. Wade 
et al 2000a) the LSD Stokes Q or U profile was integrated 
to compute the normalized equivalent width of the line po- 
larization using: 

Q = jQ(v)dv 

I J[l-I(v)]dv' 1 ' 

It was found that both net linear polarisation and longi- 
tudinal field measurements were sensitive to the integration 
range used to calculate Eqs. (1) and (2). Integration ranges 
were carefully chosen to include the entire line profile, while 
avoiding including excess continuum outside of the profile 
(which contributes only noise). This was accomplished ini- 
tially by selecting limits based on the apparent extent of 
the wings of the Stokes / profile. We then evaluated visually 
if any significant polarised flux was located outside of the 
limits, adjusting the integration limits as necessary. 

As was the case with Stokes V, LSD Stokes I profiles 
associated with Stokes Q and U profiles were locally re- 
normalised to a continuum level of 1.0 before evaluation of 
Eq. (2). 

Longitudinal field measurements are phased according 
to the ephemerides given in Table. 1. We verified that all 
adopted periods were sufficiently precise that no significant 
relative phase uncertainties exist. A harmonic curve was fit- 
ted by least-squares to the phased longitudinal field data. 
The degree of the harmonic function which yielded the low- 
est reduced x 2 , while still providing a significant and sys- 
tematic improvement to the fit, was chose as the "best" fit. 
The results of the fit are shown in Table [2] It should be 
noted that whilst the fits to the longitudinal field variations 
help compare different data sets, quantify the dispersion of 
the data and the verify the accuracy of the error bars, their 
parameters and degree do not necessarily have any easily- 
interpretable physical meaning. 

To allow a consistent comparison between the MuSi- 
CoS and ESPaDOnS/NARVAL data, the normalisation of 



the ESPaDOnS/NARVAL data had to be performed care- 
fully. Care also had to be taken with the integration limits 
when determining both the longitudinal field and net-linear 
polarization. 

When we compared the new longitudinal field measure- 
ments with those computed from MuSiCoS spectra, imme- 
diate agreement was found between the new measurements 
and those of Wade et al. (2000b), except in the case of HD 
71866 and a 2 CVn which showed a slight discrepancy. By 
re-computing the Stokes V MuSiCoS LSD profiles for these 
stars with the new masks, the observations were brought into 
agreement. It is possible that if both HD 71866 and a 2 CVn 
have sufficiently peculiar abundances, a general Ap star line 
mask (with arbitrarily enhanced and depleted abundances 
of particular elements, as used in the original analysis of 
Wade et al. 2000b) may result in an underestimation of the 
longitudinal field. Because each mask used in this study is 
"tailor-made" to the chemical abundances of the respective 
star, we are less susceptible to this problem and thus mea- 
sure slightly increased longitudinal field values. 

An example of the comparison between the new mea- 
surements, and those corresponding to the the two sets of 
MuSiCoS LSD profiles, are shown in Section 6.6 (Fig. 27) for 
a 2 CVn. This figure demonstrates that overall the agreement 
between the new measurements and the re-computed Mu- 
SiCoS measurements are acceptable. Using the new masks 
for the other stars makes little difference to the longitudi- 
nal field measurements, so the original measurements were 
kept for comparison. In addition good agreement can be 
seen between LSD profiles obtained with MuSiCoS and ES- 
PaDOnS/NARVAL observations. Fig. [T0l shows that LSD 
profiles both from MuSiCoS and the current observations 
are consistent for an identical phase of observation in the 
case of HD 32633 and HD 112413, with the new data show- 
ing much reduced noise levels. 

It was quickly seen that both the ESPaDOnS and NAR- 
VAL data were consistent with one another. An excellent ex- 
ample of this can be seen in Fig. [9] where two observations 
of HD 32633 are shown, one at a phase of 0.689 (taken with 
NARVAL) and one with a phase of 0.700 (taken with ES- 
PaDOnS). The detailed agreement between the Stokes pro- 
files indicates that both instruments appear to be perform- 
ing extremely consistently with one another. Another test 
of consistency is described in Section 6.2 and 6.6 , where for 
certain phases, the ESPaDOnS/NARVAL data for a 2 CVn 
and HD 32633 have been convolved to the same approximate 
resolution as MuSiCoS (R=35000) and then compared to 
MuSiCoS observations of a similar phase. Very good agree- 
ment can be seen in individual lines. 

As proposed by Landi Degl'Innocenti et al. (1982), as- 
sumed by Landolfi et al, (1993) and Leroy et al. (1993, 
1995), and finally confirmed by Wade et al. (2000b), the 
phenomenon of broadband or net linear polarisation results 
from differential saturation of n and a Zeeman components 
of individual spectral lines. Broadband polarisation mea- 
surements, such as those reported by Leroy (1995), effec- 
tively average the signal from the net polarised spectral lines 
with regions of continuum (which may be unpolarised, or 
polarised due to e.g. interstellar polarisation). In contrast, 
net polarisation measurements obtained using LSD measure 
only the net polarisation in lines. Therefore the measure- 
ments are not equivalent. In particular, Wade et al. (2000b) 
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Figure 7. Example abundance fit for HD 32633, with a synthetic spectrum calculation of solar abundance show with a dashed line, and 
the spectrum computed with the inferred abundances used in the LSD mask shown with a dot-dashed line. 
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Figure 9. Two observations of HD 32633 obtained at close 
phases, one with NARVAL (Phase 0.689, taken 17 Dec 2006) 
and one with ESPaDOns (Phase 0.700, taken 26 Jan 2008). The 
signal-to-noise ratio is ss 400 and ss 700 respectively. Excellent 
agreement between the signatures can been seen, illustrating the 
consistency between the two instruments. 
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Figure 10. Comparison between LSD profiles from ES- 
PaDOnS/NARVAL observations and MuSiCoS observations 
taken at near identical phases for HD 32633 and HD 112413 (Mu- 
SiCoS data dashed lines, ESPaDOnS/NARVAL solid lines). Good 
agreement can be seen. 



found that it was necessary to arbitrarily scale and shift 
the LSD measurements relative to the broadband measure- 
ments in order to bring them into agreement. Here we adopt 
a similar procedure. 
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Table [2] also shows the reduced x 2 of longitudinal field 
measurements taken from the null spectra. A value close to 
1 indicates that the null spectra and their error bars are con- 
sistent with zero longitudinal field and that the instrument 
is performing as expected. We see for HD 4778, HD 32633, 
49 Cam and HD71866, this is not the case. These high val- 
ues are a result of a handful of non-zero null measurements 
obtained during the December 2006 observations with NAR- 
VAL and very small number with obtained ESPaDonS, in 
which the null LSD profiles contained weak signatures (as 
discussed in Section 3). 

In addition, for these stars, we find that the error bars 
on the Stokes V longitudinal field measurements are larger 
than the null longitudinal field error bars, typically by a 
factor of 3-4. This is significant, because it tells us that we 
begin to see that the Stokes V longitudinal field measure- 
ments are no longer limited by photon noise. Rather, the 
larger error bars for the longitudinal field in Stokes V are 
telling us that the LSD model is failing to fit the V spec- 
trum within the Libre-ESpRIT computed error bars, likely 
due to blends and the limitation of the weak-field approx- 
imation. This only happens when a significant number of 
individual lines show very strong Stokes V signatures, and 
(as demonstrated by our ability to fit the longitudinal field 
curves within the error bars) is fully compensated for in the 
error bar calculation for Stokes V. As discussed by Wade 
et al. (2000b) this compensation is achieved by scaling the 
photon noise statistical error bars of LSD profiles by a fac- 
tor equal to the square root of the reduced % 2 ■ Wade et al. 
(2000b) describes that this scaling is almost always required 
for Stokes V in Ap stars. If we recompute the reduced x 2 
calculation again for the null spectra, but use the Stokes 
V error bars instead of the null error bars, the resulting 
reduced % 2 is consistent with a zero longitudinal field (for 
example in the case of HD 32633, the reduced \ 2 goes from 
5.42 to 0.36). 



6 RESULTS FOR INDIVIDUAL STARS 

In the following section, the results for each star will be 
discussed, including figures showing both the Stokes profile 
variation in selected individual lines and in the LSD pro- 
files. Also the variation of the longitudinal field and in some 
cases the variation in net linear polarisation as a function of 
phase will be shown for each target. Results are compared 
with those of Wade et al. (2000a) and Leroy (1995) where 
possible. In the cases where the net linear polarisation was 
close to or consistent with zero at all phases (such as HD 
32633 and HD 40312), these plots have been omitted. 



6.1 HD 4778 

HD 4778 is the only star in our sample that was not studied 
by Wade et al. (2000a & 2000b). It was classified by Renson 
and Manfroid (2009) as an AlCrSrEu star. We obtained 
Stokes V, Q and U spectra for 7 rotational phases (although 
no Stokes Q observation was obtained at phase 0.366). Using 
the method described in Sect. 5, the projected rotational 
velocity (usini) was determined to be 36 ±2 kms -1 which 
is a little larger than the Abt & Morrell (1995) value of 
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Figure 12. Longitudinal field measurements for HD 4778 ob- 
tained with ESPaDOnS/NARVAL (shown by filled circles), com- 
pared with those obtained by Bohlender (1989). A 2nd order 
Fourier fit to the ESPaDOnS/NARVAL data is given by the solid 



33 kms 1 . The data have been phased according to the 
ephemeris reported by Leone et al. (2000): 



JD = 2446674.006 + 2 d 56171 • E. 



(3) 



The individual spectral lines show strong signatures in 
Stokes V, and although variation in Stokes Q and U can be 
seen in individual lines in this star (an example is shown for 
the Fe n A 5018 in Fig. Ill[) the signatures in Stokes Q and 
U are relatively weak. Signatures are more prominent in the 
LSD profiles shown in Fig. 1111 Stokes / shows variability in 
the Fe n 4923, 5169 and 5018 A lines (as shown in Fig. fTTj) . 

Both Be and net linear polarisation measurements have 
been obtained for each phase and are reported in Table [4] 
The average uncertainly on the longitudinal field measure- 
ment is 20 gauss. The LSD profiles are illustrated in Fig. 1111 
and the longitudinal field curve shown in Fig. 1121 The 2nd 
order Fourier fit to this curve, using only the new measure- 
ments gives a reduced % 2 of 1.41. Also included in Fig. ll2l are 
measurements reported by Bolender (1989) for comparison. 
The slight apparent phase shift could be due to a cumula- 
tive error in phase over an almost 30 year span caused by 
the period uncertainty. 

Leroy et al. (1995) obtained broadband linear polarisa- 
tion measurements of HD 4778; their values are compared 
with the measurements obtained in this work in Fig. 1 131 Rea- 
sonable agreement in Stokes Q and U can be seen between 
the two epochs of data. The ESPaDonS/NARVAL data had 
the mean subtracted and then were scaled by 0.08 to bring 
them into agreement with those of Leroy et al. (1995). 



6.2 HD 32633 

HD 32633 is a fairly broad-lined B9p star with a strong, non- 
sinusoidal longitudinal magnetic field variation. We obtained 
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Figure 11. Top: Variation of Stokes I, Q, U and V LSD profiles for HD 4778. Rotational phase increases from top to bottom, represented 
on the vertical axis. A scaling factor of 25 was used for both the Stokes Q and U observations and a factor of 5 was used for Stokes V. 
Bottom: Variation of Stokes /, Q, U and V profiles for HD 4778 in the Fe II A 5018 line (scaled by a factor of 20 for Stokes Q and U and 
a factor of 4 for Stokes V). 
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Stokes V, Q and U profiles for 20 rotational phases (Stokes 
Q and U observations are missing at phase 0.086). Using 
the method described in Sect. 5, the projected rotational 
velocity (v sin i) was determined to be 19 ± 2 kms -1 which 
agrees within uncertainty with the value adopted by Wade 
et al. (2000b). 

The data have been phased according to the ephemeris 
of Adelman (1997b): 

JD = 2437635.200 + 6 d 43000 • E. (4) 

Very strong and variable Stokes Q and U polarisation 
signatures can be seen in individual lines in the spectra of 
this star, an example of which is shown for the Fe n 5018 
A line in Fig. 1151 Interestingly, the shape of the Stokes Q 
and U profiles change relatively little as a function of phase 
(in individual lines, as well as LSD profiles of Fig. I15|) . This 
could be an indication that we are "seeing" a limited part 
of the transverse geometry of the field due to the overall 
geometry of the star. It should be noted in Fig. 15 that the 
shapes of the Stokes Q and U profiles of individual Fe lines 
look quite different from the LSD profiles. This may be a 
result of differences between the Stokes profiles of weak ver- 
sus strong lines, or possibly differences in the Stokes profiles 
of lines of different chemical elements. As with all the stars 
with strong linear polarisation signatures in this study, HD 
32633 has many more individual lines showing signatures 
in the new observations than in the MuSICoS observations. 
Also Stokes / shows variability in lines such as in the Cr n 
line 4588 A and Fe n4824 and 5018 A lines (as shown in Fig. 

EJ). 

Both Be and net linear polarisation measurements have 
been obtained for each phase and are reported in Table 2] 
The longitudinal field values obtained by Borra & Land- 
street (1980) and Wade et al. (2000b) are compared to the 
new measurements in Fig. 1141 Very good agreement can be 
seen between the three epochs of data, with the clear re- 
duction in uncertainties shown by the smaller error bars as- 
sociated with the new data. The reduced \ 2 OI the Fourier 
fit to the field curve is 1.95, with an order of fit of 3. The 
average uncertainty of the longitudinal field measurements 



Figure 14. Longitudinal field measurements for HD 32633 ob- 
tained with ESPaDOnS/NARVAL (shown by filled circles), com- 
pared with those obtained by Wade et al. (2000b) with MuSi- 
CoS (shown by filled diamonds) and measurements by Borra & 
Landstreet (1980) (shown as triangles). A good agreement can 
be seen, confirming consistency between the instruments and the 
improvement in data quality is evidenced by the smaller error 
bars associated with the ESPaDOnS/NARVAL measurements. A 
3rd order Fourier fit to the ESPaDOnS/NARVAL data is given 
by the solid curve. 



is 26 gauss. The netlinear polarisation plot is not shown for 
HD 32633, due to the fact the measurements show little or 
no variation and are consistently null at all phases. In addi- 
tion for certain phases, the ESPaDOnS/NARVAL data for 
HD 32633 have been convolved to the same approximate 
resolution as MuSiCoS (R=35000) and then compared to 
MuSiCoS observations of a similar phase. Within the lim- 
its imposed by noise, very good agreement can be seen in 
individual lines between the two data sets in Fig. 1161 

6.3 HD 40312 - 6 Aur 

9 Aur is a broader-lined AOp star with a weak magnetic 
field. We obtained Stokes V, Q and U profiles for 7 rota- 
tional phases. Using the method as described above, the 
projected rotational velocity (vsini) was determined to be 
53 ± 1 kms -1 which agrees within uncertainty with the 
value adopted by Wade et al. (2000b) 

All magnetic measurements have been phased according 
to the ephemeris of Wade et al. (2000a) : 

JD = 2450001.881 + 3 d 61860 ■ E. (5) 

The individual line profiles show weak signatures in 
Stokes V; marginal Stokes Q and U signatures are visible 
at phase 0.7-0.8 (an example for the Fe II 5018 A line is 
shown in Fig. I17[) . Signatures are slightly more prominent 
in the LSD profiles. Stokes I exhibits small variations in the 
Fe II 5018 A and very subtle variation in the Fe II line 4923 
ABecause of the weak Stokes Q and U signatures in individ- 
ual lines in the spectra of this star, it would be a challenging 
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Figure 15. Top: Variation of Stokes /, Q, U and V LSD profiles for HD 32633. Rotational phase is represented on the vertical axis. With 
a scaling factor of 20 was used both the Stokes Q and U observations. Bottom: Variation of Stokes I, Q, U and V profiles for HD 32633 
in the 5018 Fe II line (scaled by a factor of 10 for Stokes Q and U and 2 for Stokes V). 



MDI of Ap Stars I. Observations 15 



HD 40312 Stokes 



-100 -50 50 100 150 
Velocity (km/s) 



HD 40312 Stokes Q 



-100 -50 50 100 150 
Velocity (km/s) 



Stokes Q 





0.01 7 



It A 



0.296 





i?F : 



0.587 



0.710 



0.789 




HD 40312 Stokes U 



0.017 



0.222 



0.296 



0.71 
0.789 
j 0.924 



HD 40312 Stokes V 



-100 -50 50 100 150 
Velocity (km/s) 



-100 -50 50 100 150 
Velocity (km/s) 







5017 5018 5019 5020 



5017 5018 5019 5020 



5017 5018 5019 5020 



5017 5018 5019 5020 



Figure 17. Top: Variation of Stokes I,Q,U and V LSD profiles for HD 40312. Rotational phase is represented on the vertical axis. A 
scaling factor of 200 was used for both the Stokes Q and U observations, while a factor of 50 was used for Stokes V. Bottom: Variation 
of Stokes 7, Q, U and V profiles for HD 40312 in the Fe II 5018 line (scaled by a factor of 80 for Stokes Q and U and 40 for Stokes V). 



16 J.Silvester et at 



Phase 0.393 Phase 0.550 




5016 5018 5020 5022 5016 5018 5020 5022 

Wavelength ( A ) Wavelength ( A ) 



Figure 16. Comparison between two phases of the ES- 
PaDOnS/NARVAL observations convolved to a resolution of 
35000 and the MuSiCoS data at a similar phase for HD 32633 
(MuSiCoS data dashed lines, ESPaDOnS/NARVAL solid lines). 
Good overall agreement can be seen between the two epochs of 
data. 



candidate for IQUV mapping, but would be very suitable 
for IV mapping. 

Both Be and net linear polarisation measurements have 
been obtained for each phase and are reported in Table 2] 
The longitudinal field values obtained by Wade et al. (2000b) 
are compared to this work in Fig. 1181 Agreement can be seen 
between the two epochs of data. The average uncertainty 
on the new longitudinal field measurement is 14 gauss. A 
Fourier fit to the longitudinal field curve yields a reduced 
X 2 of 1.53 with an order of fit of 2. Because of the limited 
phase coverage and the small net linear polarisation mea- 
surements, a plot of net linear polarisation is not shown for 
HD 40312. 

6.4 HD 62140 - 49 Cam 

49 Cam is a fairly broad-lined FOp star with a moderately 
strong magnetic field. We obtained Stokes V , Q and U pro- 
files for 19 rotational phases. The projected rotational ve- 
locity (vsini) was determined to be 24 ± 2 kms" 1 which 
agrees within uncertainty with the value adopted by Wade 
et al. (2000b). 

All measurements have been phased according to the 
ephemeris of Adelman (1997a): 

JD = 2441257.300 + 4 d 28679 ■ E. (6) 

Variable signatures can be seen in Stokes Q and U in 
individual lines in this star, an example of which is shown 
for the Fe II 5018 A line in Fig. [T3J The signatures are even 
clearer in the LSD profiles (Fig. I19[) . The amplitudes of the 
Stokes Q and U profiles are small compared to HD 32633 
(about 50 % smaller), in both the individual lines and LSD 
profiles. In Stokes I, 49 Cam shows subtle variability in most 
lines, but as with all the stars with strong linear polarisation 
signatures in this study, many individual lines show signa- 
tures, as shown in Fig. [8] 

Both Bi and net linear polarisation measurements have 
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Figure 18. Longitudinal field measurements for HD 40312 ob- 
tained with ESPaDOnS/NARVAL (shown by filled circles), com- 
pared with those obtained by Wade et al. (2000b) with MuSiCoS 
(shown by filled diamonds). A 2nd order Fourier fit to the ES- 
PaDOnS/NARVAL data is given by the solid curve. 



been obtained for each phase and are reported in Table [4] 
The longitudinal field values obtained by Wade et al. (2000b) 
are compared to those derived in this work in Fig. 1201 Very 
good agreement can be seen between the two epochs of data. 
The average uncertainty on the new longitudinal field mea- 
surement is 14 gauss, and the Fourier fit to the new mea- 
surements gives a reduced \ 2 °f 1-29 with an order of fit of 
3. To further illustrate the quality of data, Fig. l21l shows the 
longitudinal measurements obtained from the null spectrum. 
In the case of data free from any spurious signals (e.g caused 
by instrumental polarisation effects), the longitudinal field 
in the null spectrum should be consistent with zero and this 
is clearly the case with HD 62140 and is representative of 
most observations in this sample. 

The net linear polarisation as a function of phase is 
shown in Fig [22] For both Stokes Q and U measurements 
there is reasonable agreement at most phases with the mea- 
surements reported by Wade et al. (2000b) and in most cases 
very good agreement with values reported by Leroy et al. 
(1995). The net linear polarisation was scaled by the same 
factor (0.28), double the scaling factor employed by Wade 
et al. (2000b). In addition the mean was subtracted from all 
measurements, as prescribed by Wade et al. (2000b). 

6.5 HD 71866 

HD 71866 is a sharp-lined A2p star with a moderately strong 
magnetic field. We have obtained Stokes V , Q and U profiles 
for 14 rotational phases. The projected rotational velocity 
(vsini) we determine is 15 ±2 kms -1 which agrees within 
uncertainty with the value adopted by Wade et al. (2000b). 

All measurements have been phased according to the 
ephemeris of Bagnulo et al. (1995): 

JD = 2438297.5 + 6 d 80022 ■ E. (7) 
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Figure 19. Top: Variation of Stokes I,Q, U and V LSD profiles for HD 62140. Rotational phase represented on the vertical axis. With 
a scaling factor of 12 used both the Stokes Q and U observations and a scaling factor of 2 for Stokes V. Bottom: Variation of Stokes 
I, Q, U and V profiles for HD 62140 in the 5018 Fc II line (scaled by a factor of 12 for Stokes Q and U and 4 for Stokes V). 
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Figure 23. Top: Variation of Stokes I,Q, U and V LSD profiles for HD 7f866. Rotational phase represented on the vertical axis. With 
a scaling factor of 15 used both the Stokes Q and U observations and a scaling factor of 1 for Stokes V. Bottom: Variation of Stokes 
I, Q, U and V profiles for HD 71866 in the 5018 Fc II line (scaled by a factor of 16 for Stokes Q, 16 for Stokes U and 4 for Stokes V). 
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Figure 20. Longitudinal field measurements for HD 62140 ob- 
tained with ESPaDOnS/NARVAL (shown by filled circles), com- 
pared with those obtained by Wade et al. (2000b) with MuSi- 
CoS (shown by filled diamonds). A good agreement can be seen, 
confirming consistency between the instruments and the improve- 
ment in data quality is evidenced by the smaller error bars associ- 
ated with the ESPaDOnS/NARVAL measurements. A 3rd order 
Fourier fit to the ESPaDOnS/NARVAL data is given by the solid 
curve. 



Clear signatures and strong variability can be seen in 
Stokes Q and U in the individual lines of this star, an ex- 
ample of which is shown for the Fe II 5018 A line in Fig. 
1231 The amplitudes of the Stokes Q and U profiles are large 
compared to many of the other targets in this study (in both 
the individual line and LSD profiles, see Fig. I23[) . Stokes I 
appears to vary very slightly with phase, with only small 
changes in line shape in the Fe II lines at 4923, 5018 and 
5169 A . 

Both Bi and net linear polarisation measurements have 
been obtained for each phase and are reported in Table [4] 
To get a good agreement between the two epochs of data 
required a re-analysis of the Wade et al. (2000b) data, using 
the updated and abundance specific line mask (as used for 
the ESPaDOnS/NARVAL data). It is likely that this star is 
in a temperature regime where measurements are very sen- 
sitive to the line mask chosen or because of peculiar abun- 
dances. The longitudinal field variation is illustrated in Fig. 
1241 Good agreement can be seen between the two epochs of 
data. The average uncertainty of the longitudinal field mea- 
sured from the ESPaDOnS and NARVAL spectra was 27 
gauss. A Fourier fit to the longitudinal field curve with a x 2 
of 1.02, with order of fit of 4. 

The net linear polarisation as a function of phase is 
shown in Fig 1251 The Stokes Q and U measurements ob- 
tained in this study have been scaled by a factor of 0.07 to 
bring them into agreement with the measurements of Leroy 
et al. (1995) and this scaling is consistent with that of Wade 
et al. (2000b). There is an approximate agreement at most 
phases. 

HD 71866 is an ideal candidate for Magnetic Doppler 



Figure 21. Longitudinal field measurements for HD 62140 ob- 
tained with ESPaDOnS/NARVAL for the null spectrum. 
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Figure 22. Net linear polarisation measurements (Stokes Q and 
U) for 49 Cam obtained with ESPaDOnS/NARVAL (shown by 
filled circles, scaled by a factor of 0.28), compared with those ob- 
tained by Wade et al. (2000b) with MuSiCoS (shown by filled 
diamonds) and those obtained by Leroy (1995) (shown by open 
squares). A good agreement can be seen, confirming consistency 
between the instruments and the improvement in data quality 
is evidenced by the smaller error bars associated with the ES- 
PaDOnS/NARVAL measurements. 



Imaging. However, more observations will be required to 
complete phase coverage. As with all the stars with strong 
linear polarisation in this study many individual lines show 
clear signatures. 



6.6 HD 112413 



CVn 



a 2 CVn is a fairly sharp-lined AOp star with a moderately 
strong magnetic field. We obtained Stokes V, Q and U pro- 
files at 24 rotational phases. The projected rotational ve- 
locity (vsini) was determined to be 17 ± 1 kms -1 which 
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Figure 24. Longitudinal field measurements for HD 71866 ob- 
tained with ESPaDOnS/NARVAL (shown by filled circles), com- 
pared with those obtained by Wade et al. (2000b) with MuSiCoS 
(shown by filled diamonds). A good agreement can be seen (with a 
re-anaylsis of the MuSiCoS observations), confirming consistency 
between the instruments and the improvement in data quality 
is evidenced by the smaller error bars associated with the ES- 
PaDOnS/NARVAL measurements. A 4th order Fourier fit to the 
ESPaDOnS/NARVAL data is given by the solid curve. 
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Figure 25. Net linear polarisation (Stokes Q and U) measure- 
ments for HD 71866 obtained with ESPaDOnS/NARVAL (shown 
by filled circles and scaled by a factor of 0.10, compared with 
those of Leroy et al. (1995) (shown by filled diamonds). 



agrees within uncertainty with the values adopted by Wade 
et al. (2000b) and Kochukhov and Wade (2010). 

All measurements have been phased according to the 
ephemeris of Farnsworth (1932): 

JD = 2419869.720 + 5 d 46939 ■ E. (8) 

Clear signatures and strong variability can be seen in 
Stokes Q and U in individual spectral lines of this star. Pro- 
file variations of the Fe II A 5018 line and of LSD profiles are 
shown in Fig. 1261 Very strong variation in Stokes I can be 
seen in the Fe n lines at 4923, 5018, 5169 A: Both B e and 



net linear polarisation measurements have been obtained 
for each phase and are reported in Table [4] The longitu- 
dinal field values are compared to the results of Wade et 
al. (2000b) in Fig. 1271 To get a good agreement between 
the two epochs of data required a re-analysis of the Wade 
et al. (2000b) data, using the updated and abundance spe- 
cific line mask (as used for the ESPaDOnS/NARVAL data). 
It is likely that this star is in a temperature regime where 
measurements are more sensitive to the line mask chosen or 
that the abundances are sufficiently peculiar to cause a dif- 
ference. In addition the new measurements are much more 
precise. The average uncertainty on the longitudinal field 
measurement was 27 gauss, with a Fourier fit to the longitu- 
dinal field curve yielding a reduced x 2 of 0.83, with an order 
of fit of 3. 

The net linear polarisation as a function of phase is 
shown in Fig 1281 Both Stokes Q and U measurements were 
compared with the measurements of Wade et al. (2000b), 
showing very good agreement at most phases. For certain 
phases, the ESPaDOnS/NARVAL data for a 2 CVn have 
been convolved to the same approximate resolution as Mu- 
SiCoS (R=35000) and then compared to MuSiCoS as de- 
scribed for HD 32633. Again good agreement can be seen in 
individual lines between the two data sets in Fig. 1291 

6.7 HD 118022 - 78 Vir 

78 Vir is a sharp-lined Alp star with a moderately strong 
magnetic field. We obtained Stokes V, Q and U profiles for 5 
rotational phases. The projected rotational velocity (v sin i) 
was determined to be 13±1 kms' 1 which agrees within un- 
certainty with the value adopted by Wade et al. (2000b). All 
measurements have been phased according to the ephemeris 
of Preston (1969): 

JD = 2434816.90 + 3*?7220 E., (9) 

Clear signatures and strong variability can be seen in 
Stokes Q and U in the individual lines of this star, as shown 
for the Fe n A 5018 line and LSD profiles in Fig.[H]>- Stokes 
/ appears to vary only slightly between phases, with only 
small changes in the Fe n lines at 4923, 5018 and 5169 A . 

HD 118022 was studied by Khalack and Wade (2006) 
who constrained the global magnetic field of the star and 
determined the abundance distributions of titanium and 
chromium. This was performed using the magnetic charge 
distribution method (MCD) (Gerth et al. 1997; Khalack et 
al. 2001) and lower resolution MuSiCoS spectra, which lim- 
its the number of lines that could be modeled. This star is 
an excellent target for MDI, and further observations are 
warranted to supplement the currently rather sparse phase 
coverage. 

Both Bi and net linear polarisation measurements have 
been obtained for each phase and are reported in Table [4] . 
The longitudinal field values obtained by Wade et al. (2000b) 
are compared to this work in Fig. 1311 The average uncer- 
tainty on the longitudinal field measurements was 12 gauss. 
The netlinear polarisation as a function of phase is shown in 
Fig [32] and compared with the measurements of Wade et al. 
(2000b) and Leroy (1995), with good agreement achieved at 
most phases. The netlinear polarisation measurements were 
scaled by the same factor (0.10) as reported by Wade et al. 
(2000b) 
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Figure 26. Top: Variation of Stokes I,Q, U and V LSD profiles for HD fl24f3. Rotational phase is represented on the vertical axis. A 
scaling factor of 25 was used for both the Stokes Q and U observations and a scaling factor of 2 was used for Stokes V. Bottom: Variation 
of Stokes /, Q, U and V profiles of HD 112413 in the Fc II 5018 line (scaled by a factor of 30 for Stokes Q and U and 7 for Stokes V). 
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Figure 30. Top: Variation of Stokes I,Q,U and V LSD profiles for HD If 8022. Rotational phase represented on the vertical axis. With 
a scaling factor of 25 used for both the Stokes Q and U observations and no scaling on Stokes V . Bottom: Variation of Stokes /, Q, U 
and V profiles for HD ff8022 in the 50f8 Fc II line (scaled by a factor of fO for Stokes Q and U and 3 for Stokes V). 
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Figure 27. Longitudinal field measurements for HD 112413 ob- 
tained with ESPaDOnS/NARVAL (shown by filled circles), com- 
pared with those obtained by Wade et al. (2000b) with MuSiCoS 
(shown by filled diamonds). The measurements shown without 
symbols are the result of a re-anaylsis of the original Wade et al. 
(2000a) observations with the new line mask. The improvement in 
data quality is evidenced by the smaller error bars associated with 
the ESPaDOnS/NARVAL measurements. A 4th order Fourier fit 
to the ESPaDOnS/NARVAL data is given by the solid curve. 
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Figure 28. Net linear polarisation (Stokes Q and U) mea- 
surements for HD 112413 obtained with ESPaDOnS/NARVAL 
(shown by filled circles), compared with those obtained by Wade 
et al. (2000b) with MuSiCoS (shown by filled diamonds). A good 
agreement can be seen, confirming consistency between the in- 
struments and the improvement in data quality is evidenced by 
the smaller error bars associated with the ESPaDOnS/NARVAL 
measurements. 



7 DISCUSSION AND CONCLUSIONS 

The goal of this project was to obtain a new data set in 
all four Stokes parameters for a selection of well studied 
Ap stars, with the ultimate aim to map these stars using 
Magnetic Doppler Imaging. The target list contained stars 
which span a large part of the parameter space of interest, 
with sufficient signal-to-noise ratio to not only greatly im- 
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Figure 29. Comparison between two phases of the ES- 
PaDOnS/NARVAL convolved to a resolution of R=35000 and 
the MuSiCoS data at a similar phase for HD 112413 (MuSiCoS 
data dashed lines, ESPaDOnS/NARVAL solid lines). 



prove on the previous observations, but to also be suitable 
for MDI mapping. The final selection was based primarily on 
stars already identified by Wade et al. (2000a) as promising 
candidates for such study. 

Early on in the project it was clear that both ES- 
PaDOnS and NARVAL have greatly improved the level 
of detail at which Ap stars can be studied. The resulting 
dataset obtained for this study is far superior to that ob- 
tained previously with MuSiCoS, and represents some of 
the highest resolution phase-resolved observations of Ap 
stars acquired to date. This data set has more individual 
lines showing variation, much improved signal-to-noise and 
smaller error bars associated with measurements of the lon- 
gitudinal field and net linear polarisation. The new data 
have been shown to be consistent with the previous obser- 
vations of Wade et al. (2000a) and also those of Leroy at al. 
(1995), with most targets agreeing well between the different 
epochs. 

Surprisingly we found that even when the data are of 
such high signal-to-noise and when the magnetic fields are 
strong, the LSD analysis is sensitive to the normalisation 
and the measured magnetic field is rather sensitive to the 
integration ranges chosen, with variations of sometimes on 
the order of 100 gauss with very small changes of the in- 
tegration range. A key conclusion of this work is that even 
with such high-quality data, extreme care must still be taken 
with all stages of analysis to ensure consistent results at this 
level of precision. 

Although crosstalk was originally a concern, through 
a series of experiments we have shown that it is at a level 
which should not have a significant impact on the results. By 
using observations of 7 Equ, we have seen that the highest 
level of crosstalk in Stokes Q is still within the noise and 
slightly above the noise in Stokes U (around the 5 % level.). 
It should be noted that even at these levels, this effect will 
be less significant in the broader-lined stars studied here. 

Considering this, we believe the other uncertainties as- 
sociated with the analysis techniques have a greater effect: 
normalisation, blending, line masks used for LSD and the 
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Figure 31. Longitudinal field measurements for HD 118022 ob- 
tained with ESPaDOnS/NARVAL (shown by filled circles), com- 
pared with those obtained by Wade et al. (2000b) with MuSiCoS 
(shown by filled diamonds). A 1st order Fourier fit to the ES- 
PaDOnS/NARVAL data is given by the solid curve. 
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Figure 32. Net linear polarisation (Stokes Q and U) mea- 
surements for HD 118022 obtained with ESPaDOnS/NARVAL 
(shown by filled circles, scaled by 0.10 consistent with Wade 
et al. (2000b)), compared with those obtained by Wade et al. 
(2000b) with MuSiCoS (shown by filled diamonds) and Leroy et 
al. (1995) broadband linear polarisation measurements (shown by 
open squares). Good agreement can be seen, confirming consis- 
tency between the instruments. The improvement in data quality 
is evidenced by the smaller error bars associated with the ES- 
PaDOnS/NARVAL measurements. 



choice of integration ranges used for longitudinal field mea- 
surements. But as regards the final impact of the crosstalk 
on MDI mapping, this will be discussed and addressed in 
a future paper (Silvester et al. in preparation). With these 
high quality observations, we suspect that the limitations for 
mapping will in fact come not from the data (with strong 
Stokes Q and U signatures seen in many individual lines), 
but the ability to deal with line blends within the MDI code. 
An important result of this study is the confirmed sta- 



bility of the global properties of the magnetic fields of these 
Ap stars. Over multiple epochs of observations the fields 
have remained constant, with little variation in both longi- 
tudinal field and linear polarisation measurements. In some 
cases measurements separated by over a decade still agree 
with each other within the uncertainties. By comparing MDI 
maps produced from the new observations of a 2 CVn with 
those of Kochukhov and Wade (2010), we can potentially 
test for evolution of the field geometry which may occur on 
small spatial scales. 

Considering the longitudinal magnetic field measure- 
ments, agreement was found between the new measurements 
and those of Wade et al. (2000b), with the exception of HD 
71866 and a 2 CVn which showed a slight discrepancy. By 
re-reducing the MuSiCoS data with the new masks, the ob- 
servations were brought into agreement. Whilst the general 
shapes of the netlinear polarisation variations were in agree- 
ment, one needed to invoke free parameters such as scaling 
which is consistent with what was described in Wade et al. 
(2000b). 

For any study which requires observations over multi- 
ple semesters, it is imperative that the instrument is sta- 
ble and consistent throughout the campaign. Both NAR- 
VAL and ESPaDOnS proved to be very stable instruments, 
with resolution and signal-to-noise being constant over the 
4 years of data. Indeed we have also shown the two instru- 
ments are consistent with one another with close to identical 
result from similar phases. These facts demonstrate that ES- 
PaDOnS and NARVAL are both very capable instruments, 
well suited to high-resolution four Stokes measurements of 
magnetic stars over multiple year time scales. 

One of the targets (a 2 CVn) has already been mapped 
with MDI using MuSiCoS data by Kochukhov & Wade 
(2010). HD 112413 is an ideal star for determining how much 
of an improvement the new polarimetric data could give to 
MDI mapping. To quantify this improvement and to fur- 
ther confirm consistency, the new observations of a 2 CVn 
were compared with the profiles predicted by the model 
of Kochukhov & Wade (2010). As shown in Fig. ES good 
general agreement between the new observations and the 
MuSiCoS-derived model is observed. However, the new pro- 
files show more complexity than was present in the MuSiCoS 
data, which is not fully reproduced in the current model and 
would likely require a more complex magnetic field distribu- 
tion. In addition, the Stokes V profile amplitude is signifi- 
cantly underestimated by the model at a number of phases. 

The complete phase coverage of both 49 Cam (Silvester 
et al. in preparation), a 2 CVn and HD 32633 will allow the 
completion of 4 Stokes parameter MDI maps for these stars, 
doubling the number of Ap stars studied using this tech- 
nique. Out of the remaining targets HD 4778, HD 71866 
and HD 118022 would also be a worthwhile candidates for 
MDI; conversely HD 40312 has small linear polarisation sig- 
natures, relative to the noise in their spectra, making them 
less suitable for MDI analysis. The mapping of 49 Cam is 
well underway and the results will be presented in Paper II 
(Silvester et al. in preparation). 
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Figure 33. Comparison betwen the new ESPaDOnS/NARVAL data for the Fe II 5018 line (shown by black points) with the final model 
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Table 3. Log of spectropolarmetric observations, where in the 2nd to last column, an observation is denoted by a star and a missing 
observation denoted by a dash. In the final column E = ESPaDOnS and N = NARVAL. 



Object 


Date 


JD 


Phase 


iexp 


S/N ratio 


Stokes 


Instrument 




(UT) 


(2450000 +) 




(s) 


(spectral pixel -1 ) 


V QU 




HD 4778 


HD 4778 


2006 Nov 30 


4070.746 


0.423 


3200/3200/3200 


152/709/638 


-k * * 


E 




2006 Dec 04 


4074.734 


0.980 


2600/2600/2600 


281/667/858 


★ * * 


E 




2006 Dec 06 


4076.744 


0.764 


2000/2000/2000 


163/647/708 




E 




2006 Dec 19 


4089.251 


0.662 


3200/3200/3200 


129/641/552 


*** 


N 




2006 Dec 20 


4090.244 


0.050 


3200/3200/3200 


220/507/598 




N 




2006 Dec 21 


4091.260 


0.446 


3200/3200/3200 


158/585/639 


* ★ ★ 


N 




2008 Jan 09 


4475.394 


0.366 


2800/-/2800 


137/-/585 




N 


6 Aur (HD 40132) 


HD 40312 


2006 Dec 04 


4074.925 


0.587 


140/140/140 


1231/1303/1350 




E 




2006 Dec 05 


4076.147 


0.924 


52/52/52 


778/781/787 




E 




2008 Jan 07 


4473.417 


0.710 


240/240/240 


783/835/836 




N 




2008 Jan 09 


4475.268 


0.222 


240/240/240 


757/865/1056 


★ ★ ★ 


N 




2008 Jan 22 


4489.002 


0.017 


240/240/240 


1377/1225/1601 




E 




2008 Jan 23 


4490.012 


0.296 


160/160/160 


1296/906/1234 




E 




2008 Jan 25 


4491.795 


0.789 


160/160/160 


1285/1501/1468 




E 


HD 32633 


HD 32633 


2006 Nov 29 


4070.098 


0.972 


3200/3200/3200 


588/408/247 


★ * * 


E 




2006 Nov 30 


4070.830 


0.086 


3200/-/- 


646/- /- 


* 


E 




2006 Dec 04 


4074.848 


0.710 


3200/3200/3200 


789/806/776 


★ ★ ★ 


E 




2006 Dec 05 


4075.094 


0.749 


3200/3200/3200 


823/754/637 


★ ★ ★ 


E 




2006 Dec 06 


4076.085 


0.903 


3200/3200/3200 


820/755/696 




E 




2006 Dec 12 


4082.577 


0.916 


1600/1600/1600 


417/430/427 


*** 


N 




2006 Dec 13 


4083.558 


0.071 


3200/3200/3200 


379/398/426 




N 




2006 Dec 15 


4085.537 


0.379 


3200/3200/3200 


554/418/520 




N 




2006 Dec 17 


4087.528 


0.689 


3200/3200/3200 


359/414/360 


★ ★ ★ 


N 




2006 Dec 19 


4089.430 


0.984 


3200/3200/3200 


556/553/564 




N 




2006 Dec 20 


4090.419 


0.138 


3200/3200/3200 


508/455/471 




N 




2006 Dec 21 


4091.434 


0.296 


3200/3200/3200 


578/579/593 


* ★ ★ 


N 




2008 Jan 09 


4475.482 


0.017 


3200/3200/3200 


550/538/527 




N 




2008 Jan 24 


4490.736 


0.390 


2400/2400/2400 


763/788/785 




E 




2008 Jan 25 


4491.741 


0.546 


2600/2600/2600 


755/766/745 




E 




2008 Jan 26 


4492.729 


0.700 


2000/2000/2000 


672/752/726 


*** 


E 




2008 Jan 27 


4493.730 


0.855 


2000/2000/2000 


520/463/524 


* ★ ★ 


E 




2009 Jan 17 


4849.432 


0.181 


3200/3200/3200 


578/611/618 




N 




2009 Oct 10 


5115.055 


0.481 


2000/2000/2000 


659/661/695 




E 




2010 Feb 27 


5254.743 


0.209 


2000/2000/2000 


720/738/737 




E 




2010 Mar 08 


5263.774 


0.613 


2000/2000/2000 


612/600/580 




E 


49 Cam (HD 62140) 


HD 62140 


2006 Dec 10 


4080.681 


0.381 


2800/2800/2800 


685/709/735 


* * * 


N 




2006 Dec 12 


4082.653 


0.841 


3200/3200/3200 


673/757/764 


* * * 


N 




2006 Dec 13 


4083.684 


0.081 


3200/3200/3200 


518/618/599 


* * * 


N 




2006 Dec 14 


4084.595 


0.294 


3200/3200/3200 


666/505/586 




N 




2006 Dec 17 


4087.715 


0.017 


1600/1600/1600 


148/200/182 


★ * * 


N 




2006 Dec 19 


4089.557 


0.451 


3200/3200/3200 


668/628/605 




N 




2006 Dec 20 


4090.542 


0.681 


3200/3200/3200 


656/680/664 




N 




2006 Dec 21 


4091.557 


0.918 


3200/3200/3200 


774/764/743 


★ ★ ★ 


N 




2008 Jan 09 


4475.706 


0.520 


2400/2000/2000 


612/571/588 




N 




2008 Jan 24 


4491.003 


0.089 


2600/2600/2600 


957/1041/1029 




E 




2008 Jan 27 


4493.989 


0.786 


2600/2600/2600 


709/836/692 


★ ★ ★ 


E 




2008 Mar 24 


4550.724 


0.026 


2000/2000/2000 


951/929/956 




E 




2008 Mar 25 


4551.723 


0.260 


2000/2000/2000 


939/944/936 




E 




2008 Apr 24 


4581.349 


0.174 


3200/3200/3200 


605/740/655 




N 




2009 Jan 14 


4846.625 


0.056 


3200/3200/3200 


465/363/308 




N 
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Table 3. continued. 



Object 


Date 


JD 


Phase 


t-exp 


S/N ratio 


Stokes 


Instrument 




(UT) 


(2450000 +) 




(s) 


(spectral pixel - ) 


V Q U 






2009 Jan 15 


4847.620 


0.288 


3200/3200/3200 


709/723/735 


-k-k-k 


N 




2009 Sep 10 


5085.125 


0.677 


2000/2000/2000 


758/853/817 


-k-k-k 


E 




2010 Feb 23 


5250.873 


0.348 


2000/2000/2000 


860/782/809 


-k-k-k 


E 




2010 Feb 24 


5251.972 


0.604 


2000/2000/2000 


837/800/841 


■k-k-k 


E 


HD 71866 


HD 71866 


2006 Dec 15 


4085.756 


0.723 


2400/2400/2400 


527/544/534 


■k-k-k 


N 




2006 Dec 19 


4089.756 


0.312 


2000/2200/2200 


507/477/498 


-krkrk 


N 




2006 Dec 20 


4090.751 


0.457 


2600/2600/2600 


503/466/513 


-krkrk 


N 




2006 Dec 21 


4091.748 


0.603 


2800/2800/2800 


604/649/645 


■k-k-k 


N 




2007 Mar 07 


4167.957 


0.819 


2000/2000/2000 


738/748/748 


-k * * 


E 




2007 Mar 08 


4168.975 


0.969 


2000/2000/2000 


778/779/549 


*** 


E 




2007 May 06 


4227.424 


0.554 


3200/3200/3200 


549/537/518 


* * * 


N 




2008 Jan 07 


4473.530 


0.746 


3000/3000/3000 


495/423/418 


■k-k-k 


N 




2008 Jan 09 


4475.566 


0.056 


2800/2800/2800 


522/612/633 


-k-k-k 


N 




2008 Jan 26 


4492.827 


0.592 


2000/2000/2000 


699/734/651 


-k-k-k 


E 




2008 Jan 27 


4493.829 


0.740 


2000/2000/2000 


448/453/530 


-k-k-k 


E 




2008 Apr 25 


4582.411 


0.757 


3000/3000/3000 


600/444/524 


-k-k-k 


N 




2009 Jan 10 


4842.734 


0.040 


2400/2400/2400 


556/579/570 


-k-k-k 


N 




2010 Jan 28 


5224.799 


0.228 


2000/2000/2000 


685/713/708 


■k-k-k 


E 


HD 112413 


a 2 CVn 


2006 Dec 04 


4075.158 


0.977 


240/240/240 


1536/1566/1546 


-krkrk 


E 




2006 Dec 05 


4076.159 


0.160 


120/120/120 


1187/1190/1184 


■k-k-k 


E 




2006 Dec 12 


4082.760 


0.368 


240/240/240 


967/910/867 


-k-k-k 


N 




2006 Dec 19 


4089.663 


0.630 


240/240/240 


811/867/828 


-k-k-k 


N 




2006 Dec 20 


4090.647 


0.810 


240/240/240 


772/811/812 


-k-k-k 


N 




2007 Mar 02 


4163.087 


0.053 


120/120/120 


1020/1038/1049 


-k-k-k 


E 




2007 Apr 23 


4214.514 


0.457 


240/240/240 


515/572/580 


■k-k-k 


N 




2007 Apr 23 


4214.635 


0.479 


240/240/240 


528/512/503 


-k-k-k 


N 




2007 Apr 25 


4216.523 


0.824 


240/240/240 


796/781/841 


-k-k-k 


N 




2007 May 06 


4227.516 


0.834 


240/240/240 


918/764/875 


-k-k-k 


N 




2007 May 07 


4228.467 


0.007 


400/400/400 


589/494/795 


-k-k-k 


N 




2008 Jan 22 


4489.024 


0.646 


240/240/240 


1482/1555/1557 


■k-k-k 


E 




2008 Jan 23 


4490.031 


0.830 


240/240/240 


1564/1470/1569 


-k-k-k 


E 




2008 Mar 23 


4549.044 


0.621 


120/120/120 


1086/1133/1088 


-k-k-k 


E 




2008 Mar 24 


4550.850 


0.951 


120/120/120 


1011/947/1035 


-k-k-k 


E 




2008 Mar 25 


4551.891 


0.141 


120/120/120 


1190/1181/1196 


-k-k-k 


E 




2008 Apr 24 


4581.460 


0.548 


240/240/240 


931/896/906 


■k-k-k 


N 




2008 Apr 25 


4582.574 


0.752 


240/240/240 


1082/1067/1077 


-k-k-k 


N 




2008 Jul 26 


4672.734 


0.236 


120/120/120 


1010/877/1015 


•k-k-k 


E 




2008 Aug 24 


4703.724 


0.901 


120/120/120 


370/226/357 


-k-k-k 


E 




2008 Aug 24 


4703.736 


0.903 


120/120/120 


361/338/464 


-k-k-k 


E 




2009 Jan 11 


4843.728 


0.500 


240/240/240 


767/797/825 


■k-k-k 


N 




2009 Jan 13 


4846.143 


0.940 


120/120/120 


1026/1085/1039 


-k-k-k 


E 




2010 Jan 26 


5223.163 


0.873 


240/240/240 


1448/1458/1393 


•k-k-k 


E 




2010 Jun 04 


5351.852 


0.402 


240/240/240 


1665/1637/1589 


-k-k-k 


E 




2010 Jun 20 


5367.761 


0.311 


240/240/240 


1570/1557/1613 


-k-k-k 


E 




2010 Jun 22 


5369.856 


0.694 


240/240/240 


1072/983/972 


-k-k-k 


E 


78 Vir (HD 118022) 


78 Vir 


2007 Apr 23 


4214.570 


0.623 


800/800/800 


425/369/395 


•k * * 


N 




2007 Apr 25 


4216.580 


0.162 


1000/1000/1000 


327/378/346 




N 




2007 May 06 


4227.503 


0.095 


1400/1400/1400 


681/821/741 


★ ★ ★ 


N 




2008 Apr 24 


4581.520 


0.211 


1200/1200/1200 


587/634/588 


★ ★ ★ 


N 




2008 Apr 25 


4582.562 


0.490 


1200/1200/1200 


891/928/944 


★ ★ ★ 


N 
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Table 4. Longitudinal magnetic field and net linear polarisation measurements of magnetic A and B stars, measured from LSD Stokes 
V profiles and Stokes Q and U profiles. 



JD - 2450000 


Phase 


B e ± OB (G) 


Q/I ± a Q 


U/I ± a v 


HD 4778 


4070.746 


0.423 


1584 ± 22 


0.002535 ± 0.000322 


-0.002785 ± 0.000348 


4074.734 


0.980 


-1299 ±15 


-0.000206 ± 0.000341 


-0.001108 ±0.000312 


4076.744 


0.764 


-613 ±18 


-0.004080 ± 0.000390 


-0.003208 ± 0.000370 


4089.251 


0.662 


191 ±21 


0.001527 ± 0.000436 


0.007873 ± 0.000454 


4090.244 


0.050 


-1299 ± 15 


0.003215 ± 0.000516 


0.000027 ± 0.000477 


4091.260 


0.446 


1560 ± 22 


-0.003398 ± 0.000387 


0.003088 ± 0.000355 


4475.394 


0.366 


1415 ± 24 




-0.000771 ± 0.000424 


6 Aur (HD 40132) 


4074.925 


0.587 


284 ± 11 


-0.000802 ± 0.000414 


0.000139 ±0.000411 


4076.147 


0.924 


— 119 ± 17 


-0.000159 ± 0.000810 


0.001076 ± 0.000817 


4473.417 


0.710 


232 ± 17 


0.000107 ± 0.000736 


0.001391 ± 0.000738 


4475.268 


0.222 


43 ± 14 


-0.001876 ± 0.000622 


0.000512 ± 0.000534 


4489.002 


0.017 


-184 ±11 


-0.001646 ± 0.000456 


0.000054 ± 0.000389 


4490.012 


0.296 


108 ±9 


0.001905 ± 0.000388 


0.000731 ± 0.000417 


4491.795 


0.789 


156 ± 10 


-0.001068 ± 0.000446 


-0.000547 ± 0.000460 


HD 32633 


4070.098 


0.972 


-3641 ± 34 


0.000794 ± 0.000830 


-0.000148 ± 0.001156 


4070.830 


0.086 


-953 ±31 


— 


— 


4074.848 


0.710 


-2281 ± 16 


0.000290 ± 0.000383 


0.000450 ± 0.000384 


4075.094 


0.749 


-2835 ± 20 


0.000576 ± 0.000406 


0.000601 ± 0.000457 


4076.085 


0.903 


-4387 ± 28 


0.000267 ± 0.000384 


0.000267 ± 0.000401 


4082.577 


0.916 


-4245 ± 30 


-0.000343 ± 0.000688 


-0.000122 ± 0.000696 


4083.558 


0.071 


-1174 ± 33 


0.000590 ± 0.000790 


-0.000113 ± 0.000743 


4085.537 


0.379 


423 ± 11 


-0.000766 ± 0.000632 


-0.001007 ± 0.000632 


Jfl07 COO 

4087.528 


0.689 


—2165 ± 18 


0.000202 ± 0.000788 


—0.000962 ± 0.000907 


4089.430 


0.984 


-3170 ± 35 


-0.000083 ± 0.000536 


-0.000059 ± 0.000525 


4090.419 


0.138 


352 ± 32 


0.000343 ± 0.000699 


0.000409 ± 0.000673 


4091.434 


0.296 


1388 ± 18 


-0.001486 ± 0.000549 


-0.000241 ± 0.000539 


4475.482 


0.017 


-2753 ± 46 


-0.000062 ± 0.000569 


-0.000249 ± 0.000573 


4490.736 


0.390 


404 ± 12 


0.001182 ±0.000388 


0.000729 ± 0.000395 


4491.741 


0.546 


-1366 ± 17 


0.000510 ± 0.000441 


0.000231 ± 0.000436 


4492.729 


0.700 


-2185 ±20 


0.000384 ± 0.000449 


0.000414 ± 0.000449 


4493.730 


0.855 


-4256 ± 40 


0.000263 ± 0.000642 


0.000013 ± 0.000558 


4849.432 


0.181 


898 ± 18 


0.000194 ± 0.000515 


-0.000166 ± 0.000509 


5115.055 


0.481 


-925 ± 13 


0.000590 ± 0.000477 


0.000413 ± 0.000449 


5254.743 


0.209 


1273 ± 19 


0.000285 ± 0.000413 


-0.000176 ± 0.000422 


5263.774 


0.613 


-1559 ±15 


-0.000030 ± 0.000536 


0.000529 ± 0.000547 






49 C 


am (HD 62140) 




4080.681 


0.381 


-1350 ±13 


-0.001999 ± 0.000437 


-0.000067 ± 0.000459 


4082.653 


0.841 


1100 ±14 


0.001359 ± 0.000447 


-0.004501 ± 0.000452 


4083.684 


0.081 


1086 ± 14 


0.002122 ± 0.000461 


-0.003331 ± 0.000514 


4084.595 


0.294 


-819 ±19 


-0.000318 ± 0.000594 


0.001067 ± 0.000563 


4087.715 


0.017 


1253 ± 23 


-0.000198 ±0.001432 


-0.000483 ± 0.001584 


4089.557 


0.451 


-1471 ± 12 


-0.002095 ± 0.000476 


0.000024 ±0.000511 


4090.542 


0.681 


-307 ±21 


0.005206 ± 0.000519 


0.003038 ± 0.000486 


4091.557 


0.918 


1298 ± 13 


-0.001324 ± 0.000393 


-0.000273 ± 0.000438 


4475.706 


0.520 


-1358 ± 13 


-0.002141 ±0.000512 


0.001236 ± 0.000513 


4491.003 


0.089 


1112 ± 13 


0.003333 ± 0.000323 


-0.003265 ± 0.000370 


4493.989 


0.786 


598 ± 17 


0.003778 ± 0.000473 


-0.003123 ± 0.000484 


4550.724 


0.026 


1255 ± 12 


0.000867 ± 0.000322 


-0.001192 ±0.000352 


4551.723 


0.260 


-434 ± 21 


0.001381 ± 0.000389 


0.001123 ±0.000438 


4581.349 


0.174 


370 ± 21 


0.003960 ± 0.000457 


0.002064 ± 0.000531 


4846.625 


0.056 


1164 ± 14 


0.001132 ±0.000830 


-0.002572 ± 0.000928 


4847.620 


0.288 


-749 ± 20 


0.000353 ± 0.000463 


0.001085 ± 0.000497 


5085.125 


0.677 


-478 ± 20 


0.005382 ± 0.000443 


0.002041 ± 0.000419 
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Table 4. continued. 



5250.873 


0.348 


-1125 ± 14 


-0.001510 ±0.000417 


0.000311 ±0.000436 


5251.972 


0.604 


-1024 ± 13 


0.002344 ± 0.000407 


0.000870 ± 0.000401 


HD 71866 


4085.756 


0.723 


-169 ±26 


-0.001707 ± 0.000558 


-0.004811 ± 0.000603 


4089.756 


0.312 


-629 ± 20 


-0.007183 ± 0.000630 


-0.001153 ± 0.000592 


4090.751 


0.457 


-1785 ± 19 


-0.002275 ± 0.000598 


-0.003174 ± 0.000559 


4091.748 


0.603 


-1580 ± 21 


-0.003833 ± 0.000532 


0.001778 ± 0.000484 


4167.957 


0.819 


1316 ± 24 


-0.000151 ± 0.000441 


-0.003645 ± 0.000506 


4168.975 


0.969 


2223 ±21 


-0.002215 ± 0.000431 


-0.000973 ± 0.000672 


4227.424 


0.554 


-1750 ± 19 


-0.002216 ± 0.000554 


0.000866 ± 0.000535 


4473.530 


0.746 


260 ± 26 


0.001818 ± 0.000718 


-0.005288 ± 0.000736 


4475.566 


0.056 


2290 ± 22 


0.000009 ± 0.000456 


-0.000153 ± 0.000629 


4492.827 


0.592 


-1688 ± 18 


-0.003414 ± 0.000484 


0.001299 ± 0.000474 


4493.829 


0.740 


43 ±25 


0.002172 ± 0.000653 


-0.004723 ± 0.000629 


4582.411 


0.757 


461 ± 26 


0.001491 ± 0.000686 


-0.004704 ± 0.000613 


4842.734 


0.040 


2323 ± 22 


0.000023 ± 0.000475 


-0.000549 ± 0.000659 


5224.799 


0.228 


971 ± 25 


-0.004270 ± 0.000440 


0.000347 ± 0.000471 


HD 112413 


4075.158 


0.977 


-856 ± 23 


0.001780 ± 0.000494 


0.001723 ± 0.000484 


4076.159 


0.160 


-603 ± 19 


0.002334 ± 0.000525 


0.004213 ± 0.000518 


4082.760 


0.368 


624 ±21 


0.000137 ± 0.000617 


-0.001208 ± 0.000647 


4089.663 


0.630 


231 ± 23 


0.000691 ± 0.000657 


0.001269 ± 0.000687 


4090.647 


0.810 


-785 ± 19 


0.003106 ± 0.000790 


0.002424 ± 0.000788 


4163.087 


0.053 


-795 ±21 


-0.001239 ± 0.000734 


-0.000578 ± 0.000718 


4214.514 


0.457 


807 ± 18 


0.000226 ± 0.001107 


-0.001956 ±0.001125 


4214.635 


0.479 


799 ± 19 


0.001272 ± 0.000993 


-0.001996 ± 0.000981 


4216.523 


0.824 


-830 ± 20 


0.003010 ± 0.000838 


0.002400 ± 0.000780 


4227.516 


0.834 


-849 ± 19 


0.003131 ± 0.000848 


0.001945 ± 0.000738 


4228.467 


0.007 


-818 ± 24 


-0.000271 ±0.001652 


-0.001587 ± 0.000992 


4489.024 


0.646 


129 ± 23 


-0.000916 ± 0.000367 


0.001985 ± 0.000378 


4490.031 


0.830 


-851 ± 22 


0.002872 ± 0.000452 


0.002218 ± 0.000433 


4549.044 


0.621 


304 ± 22 


-0.000453 ± 0.000489 


-0.000584 ± 0.000524 


4550.850 


0.951 


-888 ± 22 


0.002450 ± 0.000779 


0.001239 ± 0.000691 


4551.891 


0.141 


-595 ±21 


0.001459 ± 0.000530 


0.004084 ± 0.000518 


4581.460 


0.548 


657 ± 20 


-0.001276 ± 0.000634 


0.001097 ± 0.000631 


4582.574 


0.752 


-595 ±21 


0.002091 ± 0.000553 


0.003702 ± 0.000545 


4672.734 


0.236 


-193 ± 23 


0.004147 ± 0.000629 


-0.001381 ± 0.000543 


4703.724 


0.901 


-854 ±31 


0.002628 ± 0.003134 


-0.000760 ± 0.001992 


4703.736 


0.903 


-907 ± 30 


0.002653 ± 0.002106 


-0.000584 ± 0.001505 


4843.728 


0.500 


-785 ± 18 


0.001821 ± 0.000779 


-0.002499 ± 0.000754 


4846.143 


0.940 


-875 ± 22 


0.002773 ± 0.000675 


-0.001509 ± 0.000701 


5223.163 


0.873 


-884 ± 19 


0.002989 ± 0.000488 


0.000467 ± 0.000503 


5351.852 


0.402 


739 ± 16 


-0.000142 ± 0.000361 


0.002013 ± 0.000375 


5367.761 


0.311 


296 ± 23 


0.001792 ± 0.000379 


-0.001181 ± 0.000359 


5369.856 


0.694 


-225 ± 21 


0.000081 ± 0.000584 


0.003735 ± 0.000598 


78 Vir (HD 118022) 


4214.570 


0.623 


-905 ± 12 


-0.005410 ± 0.000359 


0.001493 ± 0.000380 


4216.580 


0.162 


-241 ± 10 


-0.003051 ± 0.000700 


-0.005213 ± 0.000658 


4227.503 


0.095 


-127±5 


-0.001195 ±0.000454 


-0.005993 ± 0.000528 


4581.520 


0.211 


-374 ± 11 


-0.005294 ± 0.000757 


-0.003943 ± 0.000973 


4582.562 


0.490 


-957 ± 10 


-0.004111 ±0.000320 


0.002266 ± 0.000365 



